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Abstract 

A procedure to solve few-body problems is developed which is based on an expansion over a small 
parameter. The parameter is the ratio of potential energy to kinetic energy for states having not 
small hyperspherical quantum numbers, K > Kq. Dynamic equations are reduced perturbatively to 

ON 

equations in the finite-dimension subspace with K < Kq. Contributions from states with K > Kq 

o ' 

CN ' are taken into account in a closed form, i.e. without an expansion over basis functions. Estimates 

on efficiency of the approach are presented. 
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I. INTRODUCTION 



While the nuclear force cannot be treated as a perturbation, the potential energy of, say, 
a bound state of a nucleus is in fact comparable to its kinetic energy only for a limited 
number of its components. These are components having low values of the hypermomentum 
K. For all remaining components kinetic energy exceeds potential energy which allows 
solving few-nucleon problems perturbatively. 

This approach is suggested below and is a development of that of Ref. [l|] where a per- 
turbation method has been given to solve large systems of bound-state linear equations 
pertaining to the hyperspherical-hyperradial expansion. [20] The method proved to be effi- 
cient (si, |j] . The expansion parameter was the potential-to-kinetic-energy matrix element 
ratio. However, at A>3 it is the calculation of matrix elements themselves that requires a 
massive computational effort. 

The difficulty stems from a swift rise, as K increases, of the number of hyperspherical 
harmonics having the same K. (The larger is a number of particles the swifter is the rise.) 
Selection of hyperspherical harmonics to reduce the computational effort, see [l, 5, f|, is not 
efficient for A>4 bound states and is not justified in reaction calculations. In the method 
given below the difficulty is removed. No expansions are employed when K is not small. 

Recently a considerable progress in methods for solving few-body problems has been 
achieved. However, those developments have limitations not arising in the presented method. 
In particular, the well-known Green Function Monte Carlo (GFMC) method is the method 
to calculate a bound state of a system, and it does not suit to calculate reactions. (Although 
the simplest scattering problems may be considered in its frames.) Unlike this method, the 
method presented below is suitable for calculating reactions of a general type. Besides, the 
GFMC method is not convenient in that providing separate observables it does not provide 
the wave function of a bound state as an outcome of the calculation. And in the framework 
of the method described below bound state wave functions are calculated in a rather simple 
form suitable for subsequent applications. 

Recently a way was found to extend the Faddeev-Yakubovsky A=4 calculations over 
the energy range above the four-body breakup threshold |7|. However, Yakubovsky type 
calculations require too much numerical effort. The amount of calculations is less in the 
scheme below. 
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At solving few-body problems, convergence of results for calculated observables was ac- 
celerated with the help of the effective interaction methods. Such methods were developed 
in the framework of the oscillator expansion [8| and the hyperspherical expansion [9]. In 
this approach, a true Hamiltonian is replaced with some effective Hamiltonian acting in a 
subspace of only low excitations. When, formally, the latter subspace is enlarged an effective 
Hamiltonian turns to a true one. An effective Hamiltonian is constructed from a requirement 
that its ingredients, as defined in the subspace of low excitations, reproduce some properties 
of the corresponding ingredients of a true Hamiltonian in the total space. It has been shown 
8|, |9j that this, indeed, leads to improvement of convergence for observables considered. 

Higher excitations are disregarded in such type calculations. It is clear, however, that cor- 
relation effects related to higher excitations cannot be reproduced by any state vector lying 
in the allowed subspace of only low excitations. Consider e.g. the mean value, (^ol-^l^o)) 
of such an "observable" as a true Hamiltonian. It follows from the variational principle that 
an approximate state \l/o obtained with such a method provides poorer approximation to 
the true (^ol-^l^o) value than \l/o given by the simple diagonalization of a Hamiltonian in 
the same subspace of low excitations. But even the latter (^/q\H\^/q) value is a very poor 
approximation for realistic Hamiltonians. On the contrary, the method given below provides 
an approximate state vector that is apparently close to a true state vector both as to its low 
excitation component and to its high excitation component. 

And speaking of reaction calculations in the framework of the dynamic schemes employed 
below one should also take into account that a rate of convergence is determined not only by 
the properties of a Hamiltonian but also by those of the source-term q entering the equations. 
But these properties are apparently ignored at constructing effective Hamiltonians. Unlike 
this, the method developed below provides state vectors genuinely close to the true ones 
both for bound state problems and reaction problems. 

Development of efficient microscopic methods for nuclear physics is timely now because 
of a necessity to test nuclear forces derived from the effective field theory. 

The method is given in the next section, and its implementation is considered in Sec. 3. 
In Sec. 4. estimates on its efficiency are presented. In particular, the method seems to be 
promising to test realistic nuclear forces in four-nucleon reaction problems and in scattering 
problems below three-body breakup thresholds in the 5 He, 7 Li, 7 Be, and 8 Be systems. This 
opens up basically a new field. 
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II. THE SMALL PARAMETER EXPANSION 



We consider first the bound state problem 

(H-Ei)% = 0, (1) 

where H = T + V is an A-body Hamiltonian. We split the whole space of states into the 
subspaces with K < K and K > K and we denote ^\ and ^ the components of the 
solution that lie, respectively, in these subspaces. Let us denote Pk and Qk projectors 
onto the K < K and K > K subspaces, respectively. Let be the zero approximation 
eigenvalue that arises when coupling to the K > K subspace is disregarded, 

P Ko (H - Ef) *f> = 0. (2) 

The state tyf ^ is the corresponding zero approximation eigenf unction. 

Let us also define a Green function in the K > K subspace. The operators of hyper- 
rotations commute with the kinetic energy operator. Therefore, if \I/q is a state lying in 
the K > K Q subspace then T^fg also belongs to this subspace. We shall consider also the 
"correction operators" 5H specified below that have the same propety. Then 

VQ = {T + 5H-Ef ) )m Q (3) 

is a state belonging to the K > K subspace. We also assume that T + 5H possess 
only continuum spectrum. Then for any state /q in the K > K Q subspace one has 
(T + 5H — E^)/q ^ 0. Taking also into account that the operator from (jSJ) possess a 
complete set of eigenstates with K > K Q one concludes that there exists a unique state 
corresponding to any given (fQ in ([3]). Thus we can define the corresponding Green function 

G Ko = (t + SH-E^Y'qk,. (4) 

Taking into account that Pk TQ Kq = Qk TPr = we write down Eq. (pQ) in the form 

P Ko (H-E^\ = -P Ko V^, (5) 
= -G Ko [V¥t + (U- AEi) , (6) 

where 

U = V-8H, AEi = Ei - Ef . 
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When the hyperangular momentum K is sufficiently large kinetic energies of states in the 
K > K subspace are high. Then in accordance with Eq. (j3J) is "small", so that one 
may treat the second term in the right-hand side of Eq. (jHJ) as a perturbation. Thus one 
may express the component ^ in terms of ty\ perturbatively and obtain dynamic equations 
for the latter component alone. Namely, one may write 

= -rv% (7) 

where 

T = G Ko - G Ko (U - AEi)G Ko + G Ko (U - AEi)G Ko (U - AEk)G Ko (8) 

(We mention that, say, the contribution from AE{G\ here is not of the same magnitude 
as that from Gk UGk -) The dynamic equation (TjQ) may then be put in the form of an 
equation in the subspace of low excitations only, 

P Ko (T + V - VTV - Ei)^ = 0. (9) 

The quantity — VTV represents a genuine effective interaction that reproduces the effect of 
the excluded K > K subspace. 

It is convenient to treat perturbatively not only subsequent contributions to the T oper- 
ator but also the whole —VTV interaction. This leads to an expansion over powers of both 
Gk U and Gk V. We replace Gk with AGa' and we seek for ^(A) as an expansion, 

¥j(A) = £ X m ^ m \ E(X) = J2 A m 4 w) (10) 

m>0 m>0 

setting A = 1 at the end. The zero order eigenstates and eigenvalues are given by Eq. (J2j). 
Equations for higher order corrections are of the form 

P Ko {H - Ef ] )^ m) =q^ (11) 

with a given source term. Explicit form of the equations for and \l/^ 2 ' ) is 

P Ko (H - = E^ 0) + P Ko VG Ko V^°\ (12) 

P Ko (H - *f )*f > = E^f ] + E?hf + P Ko VG Ko V^ 

-P Ko VG Ko UG Ko V^f ] . (13) 



(The term from (jSJ) with AEi — > would first appear at m — 3.) 

Obviously, solutions to Eqs. (Ilip are not unique: at given lower order corrections entering 
q(m) foe solution v[o m ) is determined up to const . It is convenient to select unique 
solutions imposing e.g. the usual perturbation theory condition 

(*ri*!(A)>H*fw o) >- (14) 

This is a normalization condition since up to a numerical factor ^'(A) are determined by 
Eq. (0Q) type equations. According to the first of Eqs. (ITU]) this is equivalent to the set of 
conditions 

(^f)|^( m )) = 0, m = l,2,.... (15) 

(Use of any other condition of Eq. (I14p type whose right-hand side includes additional terms 
tending to zero at A — ► would lead merely to a different normalization of the sum 

In order Eqs. fTTTT) be self-consistent the conditions (^f^lq^) = should be fulfilled. 
This determines the energy corrections E[ m \ ( In the present consideration complications 
that arise in cases of possible non-trivial (quasi) degeneracies of levels are disregarded.) We 
get, in particular, 

m (^ m \VG K V\^ m ) 



{^f ) \VG Ko UG Ko V\^f ) ) - {^ f^VG^V^P) 
(*fW 0) > 



Since by construction Gk is positive definite the correction E^ 1 ' is negative. (Besides, one 
has in (HZD (^ 0) \VG Ko V\¥} 1] ) = -(^ 1] \H - ^ (0) |^ (1) ).) 

Thus we have dynamic equations of the Eqs. fllUp . fill I) form which do not involve 
excitations with K > Kq at all. This is an advantage since the K < Kq subspace spans only 
a finite, hopefully not too large, number of HH. 

The complementary K > Kq component of a state sought for is given by Eq. (JTJ). This 
can be represented as 

where, in particular, 

^« = -G Ko Vmf\ (19) 
V- {2) = -G Ka V^f ] + G Ko UG Ko Vy l f°\ (20) 
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The low excitation component ty\ is obtained above in the form of a hyperspherical 
expansion. This component may be stored in this form for use in the relations of the 
Eqs. (1191) . (1201) type and also in other applications. The complementary high-excitation 
component is reconstructed as a quadrature and may be applied in such a form. This 
component does not involve any expansion. 

According to the variational principle the energy eigenvalue calculated as the average 
value of a Hamiltonian over \l/^ ' ) + + is accurate up to the order m — 3. 

Let us also comment on applying of a such type expansion for calculation of reactions. 



10 



First we shall proceed in the framework of the following well-known procedure 
Consider reactions below the three-fragment breakup threshold when only two-fragment 
channels are open. Let be a continuum spectrum state and denote N the number of open 
two-fragment channels. The ansatz 

N 

*. = # + Erf+ X > (21) 

is used where and (j^p represent the "channel" states of two possible types, are 
reaction amplitudes to be determined, while X is localized and is sought for as an expansion 
over hyperspherical harmonics. In the three-nucleon case the procedure is applicable also 
above the three-nucleon breakup threshold. In this case X describes breakup to free particles 
at large distances. The equation determining X is 

(H-E)X = q, (22) 

N 

?=-0, (1) +£^f, 

3=1 

where 4 (1),(2) = (H - E)^'^ are localized states. The state X is found from Eq. (1221) up 
to reaction amplitudes /y (cf. below), 

N 

X X,~Y,fu X .r (23) 

3=1 

and fij are obtained from iV additional linear equations. 
Proceeding as above we represent X as X 1 + X h . We have 

P Ko (T + V-VTV- E) X 1 = P Ko (q - VTq) , (24) 
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X h = -r (VX l - q) . (25) 

Here T is of the form of Eq. (jSj) with being omitted and with Ei — > E in Eq. (J3J). 
Writing as above X 1 = ^2 m X^ m \ X h = Y^ m ^ one obtains equations for X 1 ^ and 
expressions for X h(jn ^ performing an expansion over Gk - One has, in particular, 

P Ko (H - E)X l ® = P Ko q (26) 

while the corresponding equations for X 1 ^'^ 2 ' and expressions for X h ^^ 

are obtained 

from Eqs. (fT2"j) . ffTHj) and (TTTJ]) . f[2"U]) with i^ 1 -^ 2 ) being omitted and with the replacements 
Ef ] E, and #^(i)>(2) -> X^ 1 ^, -> VX'^ - g. 

One may note that Eq. f[2"2"j) has a localized solution only when the amplitudes /y equal 
to their true values. Otherwise, the solution includes an admixture of cluster components. 
On the contrary, if a finite number of terms is retained in the expansion of T over powers of 
Gk then the localized solution exists at any allowing the representation (T23]) . The reason 
is that the role of higher terms in the expansion of T increases at large distances. When 
in (T2"2"]) not coincide with their true values these higher terms are responsible for description 
of the large distance clusterization (cf. below). The same occurs in the approximate way 
to solve Eq. ( 1221) applied up to now when the equation of the form ( 1261) was employed. Of 
course, this does not pose any problems. 

In the general type reaction case the method is applicable in the framework of the ap- 

rfnn 

proach [12J, [13|, HJj in which reaction observables are obtained from states \l/ that vanish at 
large distances like bound states. The approach extensively applied for perturbation induced 
reactions and proved to be very efficient. Any strong-interaction induced reactions can also 
be treated in this way. Dynamic equations for states \l/ are of the form of Eq. ( 1221) with 
X — > \I/ where q is a given state and the energy E is complex. Because of the latter the 
solution ^ is localized and it is a proper object to be found with the help of the expansion 
over Gk - The perturbative expansion to calculate = \E^ Z + is similar to that described 
above. 

The outcome of a calculation are quantities of the form Q(E) = (^/'(E)\^/(E)) where ^' 
is a state similar to \1/ for another source term q' . Reaction observables are extracted from 
&(E) in a simple way as quadratures. When it is sufficient to calculate &{E) only up to the 
first order in Gk one need not account for ty h whose contribution in Q(E) is of the second 
order. 
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III. IMPLEMENTATION 



We are dealing with the space of Jacobi coordinates or that of Jacobi momenta and we 
denote n = 3A — 3 the space dimension. Considering matrix elements in the momentum 
representation we denote tt the n-dimensional momentum vectors and II = |7f| the hyper- 
momentum. We adopt in (JSJ) SH = V(U) where V^(II) is a subsidiary interaction. We 

(7F ^ M = IP /(2m) + W) -E ' (2?) 

(7f'|C/|7f) = (7t'\V\7t) - ^(tK - k)V(I1), (28) 

where for bound states E means E?\ Here Y Kv is an orthonormalized complete set of 
hyperspherical harmonics with a given K , and ft denotes a unit vector pointed in the direction 
of ft, 7r = 7f/II, 7r' = n' /H' . The hyperangular factor entering (1271) may be represented with 
the simple expression (e.g. 

_ J{ + — fn-2\ "-2 

E Y ^(^)^(7r) = f — J ^(tt'-tt), (29) 

where CT K {p) is the Gegenbauer polynomial. 

When performing calculations in the coordinate representation we denote £ the n- 
dimensional position vectors, p = |£|, and £ = £/p. We have 

2M p 2 



T = T ? + —^ (3°) 



«i^o^(e'-o(-^)($ + ^). (3D 

Here K 2 is the hyperangular momentum operator. In this case we choose: 

SH = V(p)-T p + E, (32) 

where V(p) is a subsidiary interaction. It is convenient to represent the corresponding Gk 
as a sum of contributions from various K values, 

G Ko = 9k- (33) 
k>k 



Then 



h 2 K(K + ii - 2) - , ' 



(C\U\0 = (?\V\0 - * (n) (£' - [V(p) -T p + E}. (35) 

The choice ( 1321) is done to facilitate Monte-Carlo calculations of matrix elements. In the 
coordinate representation, Gk that contains the total T as in Eq. ( 1271) would correspond 
to qk with hyperradial Green functions varying rapidly at not small K. 

When performing calculations it is convenient to include the factor = I hi 

the Green functions, where {6^} is a complete set of spin-isospin states (c.f. [3]) with simple 
permutational properties. 

The subsidiary interactions V(p) and V"(II) can be chosen from the requirement of fast 
convergence of observables as Kq increases when one takes into account only, say, the lowest 
order corrections. Alternatively, one can minimize the ratio of the second order correction 
to the first order correction for this purpose. 

If e.g. in the coordinate representation one writes (suppressing the spin-isospin notation) 

*? m) = £ (36) 

K<K ;u 

then the equations of Eqs. (TTTT) - (fT3l) and (I2"B"]) type are of the form 

- m {i + - ^7^) « - + = ^ 

(37) 

where K < K and K' < K . 

Their right-hand sides as well as the \|> /l ( m ) components are to be calculated with the 
Monte-Carlo method. Integrands depend on high K values only via Gegenbauer polynomials 
(|29p entering (|27|) and While these polynomials are rather quickly oscillating all other 
factors in the integrands are smooth functions of coordinates or momenta. In the case m = 1 
one may simplify a calculation taking the argument of Gegenbauer polynomials as a new 
integration variable. Integration over this variable may be done with the help of the regular 
Gauss-Gegenbauer quadratures while integration over other variables that are smooth may 
be done with the Monte-Carlo method. A suitable change of variables is described in 
Appendix. This can also be done in the case of a momentum representation calculation. At 
the same time there are indications (e.g. 16| ) that direct Monte-Carlo integration may be 
suitable even at rather large K values. 

Local components are dominating components of nuclear forces derived from the effective 
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field theory. If m > 2 correction terms are retained in a calculation a reasonable simplifica- 
tion may be to account for only those local components in these terms. 

Direct solution of Eqs. (1371) in the form they are written down is hampered by the large 
centrifugal barriers K(K + n — 2)/p 2 . In the bound-state case a practical procedure is to 
expand Xku (p) over a se ^ °f functions that reduces Eqs. (1371) to linear equations. Such linear 
equation sets of a large size may efficiently be solved with a version of the method of Ref. 
1|, i.e. using an expansion over another parameter of the K^ 2 type. For complex E values 
with positive real parts and rather small imaginary parts entering Eqs. (13*7)) in the reaction 
calculations convergence of the expansion procedure to solve Eqs. ( I3~7j) is slow. And for real 
E values entering Eqs. (1371) in the other type calculations of reactions this procedure may 
lead to unphysical singularities in reaction observables. Other efficient solution methods are 
available for this purpose. 

Integral transforms <&(E) are required at sufficiently many values of complex energies E 
to perform a satisfactory inversion [l2, k], Q]- But one need not solve Eqs. (137)) for all 
these E values. A better approach is to solve these equations for a rather small number of 
E values and to obtain $(£7) for a larger set of E values via interpolation. The transforms 
<&(E) are smooth functions and this procedure is safe and accurate. 

IV. ESTIMATES 

First let us consider numerical estimates for the 4 He system in the + state. In 3j the 
accurate a-particle binding energy pertaining to a proper subset of hyperspherical harmonics 
has been compared with energies calculated approximately as follows. Only matrix elements 
of NN force {Kv\V\K'v') such that either K or K' does not exceed some Kq were retained 
in the system of equations. This approximation is equivalent to the approximation V — > Gk 
in Eq. ([91). The approximate solution thus obtained is close to that given by the m = 1 
approximation of the method described above. It accounts for corrections to binding energy 
up to m = 3. The corresponding approximate energy values as a function of Kq along with 
the accurate value are shown in Table 4 of Ref. jjj. 

At Kq = 14 the difference between the two values equals to 0.23 MeV for the NN inter- 
action with a very strong repulsive core. In this connection one needs to remember that the 
binding energy considered is a small difference between two large quantities, potential and 
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kinetic energy, which deteriorates the accuracy. This comparison refers to the case when the 
above mentioned mean field V is set to zero. Inclusion of the mean field in the calculation 
would improve the convergence. Furthermore, it is seen from Table 4 of [3| that for another 
NN force having a softer core convergence with respect to Kq is much faster. One may 
note in this connection that nuclear forces derived from the effective field theory are much 
softer than that used in the above comparison. For such forces convergence with respect to a 
maximum K value retained in a calculation is considerably faster than for phenomenological 
local realistic forces 1^]. Naturally, for effective field theory forces one may expect faster 
convergence as to Kq as well. 

Now let us perform an estimate for the case of reactions in the same system. We consider 
the above outlined approach dealing with complex energy E, and we set E = <7r + 207. 



We take 07 = 10 MeV which is a good value for e.g. electromagnetic processes [14j. We 
consider again the case K = 14 and again we shall not include the mean field V. We 
employ the so called AV4 potential that is the central component of the realistic AV18 NN 
interaction 18]. For the estimate purposes we adopt the following model. We represent the 
subspace with K < Kq with a single hyperspherical harmonic with K = 0. We represent the 
subspace with K > Kq with a single hyperspherical space-symmetric "potential" harmonic 
with K = 16. Thus in our model we deal with two coupled differential equations that have 
the form of Eqs. (J3"7j) and that correspond to the hyperspherical expansion of the state 
\& described in the preceding section. The right-hand side source term in the first of the 
equations was set to be exp(— p/0.4fm). The form of the source terms is not very important 
for the estimate and in the second of the equations the source term was set to zero. In 
general, at large p values the state \I/ contains the 3N+N cluster components, in particular. 
When K is less or about p/y3R, R being the range of the 3N cluster, in equations of 
the type we consider matrix elements may produce coherent effects. Then it would not be 
suitable for our estimates to represent the K > Kq subspace with a single hyperspherical 
harmonics. In our case, however, a typical extension in p of the state \I/ is 10 fm, c.f. below, 
and with our Kq value these effects are not relevant. The results for the quantity of interest 
<§>(E) = (^(E)l'fy(E)) are shown in Fig. 1. The curves labeled as the first and the second 
approximation correspond, respectively, to calculations in the frameworks of m = 1 and 
m = 2 approximations as described in the preceding section. We also note that the central 
force we use for the estimate has a strong repulsive core of the height of 2.7 GeV. As above, 
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FIG. 1: Integral transform of the form factor of the type of those governing reaction amplitudes 
or representing response functions. The calculation is done for the model described in the text. 
The quantity 07 was set to be 10 MeV. Exact values and the m=l and m=2 approximations are 
compared. 



one may note that use of a nuclear force derived from the effective field theory and inclusion 
of the mean field V in a calculation would improve the rate of convergence. 

In the problems we considered there exists several hundreds hyperspherical harmonics 
with K < 14 that is acceptable. A swift rise of the number of hyperspherical harmonics 
with the same K starts at K values about 10, see Table 3 from So it would be desirable 
to keep Kq at the level of 10-12 in a calculation. 

In general, conditions for applicability of the perturbation approach we consider are 

\\G Ko Uy h \\ < ||* ft ||, \\G Ko V¥\\ < ||*'||. (38) 

A better accuracy in the first of these conditions is provided reducing U via inclusion of the 
mean field V. The second of these conditions contains the interaction in the form Qk VPk - 
It is non-diagonal in hyperspherical quantum numbers and is substantially "smaller" than 
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Qk q VQk q for this reason. For bound states and for scattering states in the most interesting 
case when KeE is smaller or about \(V — V)\ the first of these conditions, for example, may 
roughly be represented as 



where p ~ A l l 2 (r 2 ) 1 / 2 is the corresponding range either of the bound state or that of the 
component X of the continuum state above, or that of the state ^ above while (V — V) is 
the average value over or X, or \l/.[21( In the case, for example, of the 7 Li two-cluster 
ground state the rate of convergence may be increased if one diminishes p separating out 
the cluster component of the state (c.f. Eq. ( I21i) . 



The value of p in ( 1391) in the case of the ^> state is about (Imfc) 1 where {hk) 2 / (2M) = E. 
If E = aR + icF! then Imfc = 2~ 1 / 2 [(o-| + a|) 1/2 - cr R } 1/2 . The choice of 07 is discussed in Q, 



It occurs that when KeE = <jr increases at given K Q the left-hand sides in Eqs. (I3"H]) 
decrease if Gk depends on E as in Eq. ( 1271) . 

We note that the separation out of two-cluster components to diminish the p value may 
also be done when calculating the \1/ quantities. At large distances the relative motion factor 
is an outgoing (neutral or Coulomb) wave with a complex wave number having a positive 
imaginary part in this case. This ensures vanishing of these components at large distances. 

In the n+n+p continuum-state case the method may be applied also above the three- 
body breakup threshold despite the fact that the component X is not localized. The reason 
is that average potential energy for the total breakup channel decreases with p as p~ 3 . But 
in the 3-body case it may be profitable to do the whole calculation in the matrix form i.e. 
to use the HH expansion for the perturbative calculation also of the K > K contributions. 

In conclusion, an area for testing realistic nuclear forces may be substantially extended 
with the help of the presented approach. For this purpose, it is required to investigate the 
issue of Monte-Carlo computing the perturbative corrections. 



The work was partially supported by RFFI, grant 07-02-01222a, and RMES, grant NS- 
3004.2008.2. 




(39) 



see also 1191. 



14 



APPENDIX 



When one takes f ■ ( as a new integration variable one needs to define the whole set of 
integration variables in a way that the integrand remains non-singular. This can be done 
e.g. as follows. Let us express a unit vector £ = {£1, . . . , £„} in terms of another unit vector 
V, 

n 

ii = ^ ] 9ijVji 
3=1 

where g^ is an orthogonal matrix such that its first column is gn = and gij is arbitrary 
otherwise. One then has £' • £ = £V j9n9ijVj = Vi- Let us parametrize the components of r) 
as follows, 

571 = cos ip, f]j = sin tp, j = 2, . . . , n, 

where Vi are components of a unit vector -0 on a hypersphere in a n — 1-dimensional subspace. 
Taking into account that 

d£ — dr) = (sin ip) n ~ 2 dvdip 
one then may rewrite e.g. integrals of the structure (F 1 \Gk \F 2 ) as 

E ^7^ T {^) f P^dpismvT-^dvdi'FKpi) 



K>K 

,1 TStTS 1 „ o\ I" 1 

C x 2 (cos^)F 2 (pO, 



^ ^ + n-2) +y(p) 



n-2 



2m p 2 

where the components of the n-dimensional unit vector £ entering i*\ are parametrized as 
follows, 

ii = §cos</? + (^tj9ij(i')vj-?J sin p. 

The integrations over dp, d£', and dv may be performed with the Monte-Carlo method 
while the remaining integration over dip may be done with the help of regular quadratures. 
Integrals at a given p over the hypersphere which pertain to the m = 1 correction are 
transformed similarly. 
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